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Abstract 


This  numerical  investigation  paremetrically  determined  how 
material  and  geometric  properties  affect  the  buckling  load 
capability  of  8-ply  curved  composite  panels  loaded  under  pure 
shear.  The  linear  bifurcation  capabilities  of  the  STAGSC-1  finite 
element  computer  code  developed  by  the  Lockheed  Palo  Alto 
Research  Laboratory  was  employed.  Results  from  this  computer  code 
were  compared  with  other  published  analytical  results  of  similar 
configuration  to  verify  boundary  conditions. 

Representitive  material  properties  for  both  graphite-epoxy 
and  Kevlar-epoxy  were  used.  The  aspect  ratios  of  0.5,  1.0,  and 
2.0  were  achieved  by  varying  the  axial  length  of  the  cylindrical 
panel.  Both  quasi-isotropic  and  symmetric  orthotropic  fiber 
orientations  were  analysed. 

Results  indicate  that  these  curved  composite  panels,  while 
loaded  with  a  shear  load  resulting  in  pure  shear  reactions  in 
flat  plates,  do  not  react  in  pure  shear.  Once  the  interior  force 
resultants  are  established,  the  load  carrying  capability  of  the 
curved  panels  follow  classical  laminated  plate  theory. 


ANALYTICAL  SHEAR  BUCKLING  INVESTIGATION  OF 
CURVED  COMPOSITE  PANELS 

I.  Introduction 

Background 

As  we  try  to  make  our  vehicles  go  further,  higher,  and 
faster  we  enevitably  run  into  one  requirement;  make  them  lighter 
and  stronger.  Many  studies  on  future  systems  have  shown  the  need 
for  lighter  structures.  Some  examples  are  the  Transatmospheric 
Vehicle,  the  Advanced  Tactical  Fighter,  and  almost  all  space 
systems.  This  is  especially  true  for  any  of  the  space  systems 
since  having  to  get  out  of  the  gravity  well  is  expensive. 

One  possible  solution  to  this  problem  is  the  use  of 
composites.  Composites  provide  a  light,  strong  material  that  can 
be  tailored  to  the  purpose  by  orienting  the  fibers  to  afford 
strength  in  desired  directions.  The  composites  also  promise 
better  fatique  life,  damage  tolarance  [10]  and  can  be  formed  more 
easily  to  complex  curves.  The  present  problem  with  composites  is 
the  fact  that  they  do  not  react  like  isotropic  materials.  The 
fact  that  their  stiffness  can  be  oriented  to  preferred  directions 
make  them  complex  to  analyes. 

As  Whitney  states  in  his  latest  paper  [13],  "Numerous  papers 
concerning  the  instability  of  laminated,  anisotropic  plates  and 
shells  can  be  found  in  the  open  literature.  Buckling  of  curved 
panels  has,  however,  received  little  attention"  This  is  even  more 
the  case  with  shear  buckling.  The  work  done  in  Refs.  [2] ,  [3] , 
and  [4]  is  concerned  with  compressive  buckling  of  curved  panels. 


Reference  [1]  deals  with  shear  buckling  of  an  entire  cylinder, 
not  a  curved  panel.  In  this  thesis,  a  parametric  study  is 
performed  on  a  curved  composite  panel  and  how  the  parameters, 
described  next,  effect  shear  buckling  loads  using  a  finite 
element  code. 

Problem  Statement 

This  study  examines  the  buckling  characteristics  of  a  simply 
supported  cylindrical  composite  panel  with  a  12  inch  radius  of 
curvature  loaded  in  pure  shear.  The  parameters  which  were  varied 
were  aspect  ratio,  fibre  orientations,  and  material  properties. 
The  values  for  these  parameters  are  presented  in  chapter  3.  The 
results  are  analysed  and  explanations  for  the  trends  are 


The  main  difference  between  bifurcation  buckling  and 
collapse  load  (as  the  terms  are  used  here)  is  that  bifurcation 
buckling  uses  the  original  stiffness  matrix  computed  for  the 
structure  through  buckling.  A  collapse  load  analysis  is  a 
nonlinear  method  that  continually  updates  the  stiffness  matrix  as 
the  structure  deforms  and  takes  on  new  shapes.  Both  techniques 
increment  the  applied  load  until  the  point  that  produces  a 
structural  instabillity  is  reached.  This  is  characterized  for  a 
linear  prebuckling  solution  when  the  structure's  equilibrium 
equations  have  multiple  solutions. 

Comparisons  of  the  results  of  bifurcation  and  collapse  load 
analysis  [8]  have  indicated  when  each  should  be  used.  The 
bifurcation  analysis  normally  predicts  a  lower  failure  load  than 
does  the  nonlinear  collapse  load  technique.  However,  in 
structures  not  containing  geometric  discontinuities,  such  as 
holes,  the  bifurcation  analysis  provides  a  fair  approximation  of 
the  collapse  load  results.  When  structures  contain  geometric 
discontinuities  the  nonlinear  results  indicate  that  the 
deformations  along  the  free  edges  of  the  hole  permit  a 
redistribution  of  the  stresses  away  from  this  free  edge 
throughout  the  structure,  thus  resulting  in  a  higher  load  bearing 
capability. 

Since  the  structures  in  this  thesis  are  simply  supported 
shells  with  no  cutouts,  the  linear  bifurcation  portion  of  the 


STAGSC-1  [6]  finite  ellement  code  was  employed.  The  theory  for 
this  portion  of  the  code  is  discussed  in  the  STAGSC-1  Theory 
section. 


Classical  Laminated  Plate  Theory 

An  understanding  of  classical  laminated  plate  theory  (CLPT) 
is  required  whenever  dealing  with  composite  structures.  This 
section  is  intended  to  provide  a  brief  overview  of  basic 
principles.  A  more  complete  treatment  can  be  found  in  any 
introductory  composite  text  such  as  [10].  CLPT  is  presented  here 
for  completness. 

We  begin  by  discussing  the  stress-strain  relations  for  a 
single  lamina  in  plane  stress.  These  relations  are  given  by 


°U  °L2  0 
=  %1  %2  ° 

0  0  %6 


where: 


(1) 


QL1=EL/(l-v12v2i) 

°L2  =V12  *2  /  {1~V12  V21  ] 

(2) 


%2  =E2/ll-^2VZl) 


Note  that  the  1  and  2  subscripts  denote  the  principal  material 
direction  (see  Fig.  1) . 

For  the  more  general  case  where  the  principal  material 
directions  do  not  align  with  the  natural  body  axes,  the  reduced 
stiffness  matrix  [Q(J]  must  be  replaced  by  the  transformed  reduced 
stiffness  marix  [Q(J]  where: 

—  4  2  2  4 

Q-j^Q-^cos  »  +  2  (Q12  +  2Q^^)  sin  0cos  0+Q22sin  ® 

_  2  2  4  4 

Ql2=  (Qh+Q22~4®66^  S^n  9cos  ®+Q^2^s^n  ®+cos 

—  4  2  2  4 

^22=^llsin  9+2  ^12+2^66^ Sin  flcos  ®+®22cos  9 

(3) 

Ql6=  (Qii"Qi2"2q66)  sin9cos^9*  (Q12"Q22+2Q66)  sin^cos* 

3  3 

?26=  (qH-q12_2Q66)  sin  ®cos*+  (q12-q22+2^66)  sin*cos  9 

®66=  (Q11+Q  -2Q12-2Q66)  sin2®cos29+Q66(sin^a+cos^d) 

Note  that  the  [QtJ]  matrix  is  fully  populated  and  the  stresses  and 
strains  in  Eq.  (1)  are  measured  along  the  natural  body  axes,  x 
and  y. 

To  expand  these  arguments  to  a  general  layered  laminate,  two 
assumptions  must  be  made.  The  first  is  that  the  bonds  between  the 
lamina  are  perfect.  That  is,  there  is  no  slipping  between 
laminae.  Secondly,  the  Kirchhof f-Love  hypothesis  which  states 


that  normals  to  the  midsurface  remain  normal  after  deformation  is 
enforced  throughout  the  laminate.  These  assumptions  lead  to  the 
strain  equations 

K 

X 

*y  '  ,4) 

«::y  , 

where  the  o  superscripts  represent  midsurface  values  and  the  k's 
are  midsurface  curvatures.  The  distance  from  the  midsurface 
perpendicular  to  the  x-y  plane  is  given  by  z  as  shown  in  Fig.  2. 

The  above  strains  and  curvatures  are  given  by  Sanders' 
kinematic  relations  [9]  for  a  cylindrical  panel  with  moderately 
large  displacements  and  rotations  of  tangents  to  the  midsurface 
by 

«  °  =u  +0.5*2+0.5*2 

A  9  A  .A. 

«  =v  „  +  (w/r)  +0 . 5*2-0 . 50 


where  u,  v,  and  w  are  the  axial,  circumferential,  and  radial 
displacements,  respectively,  of  the  panel's  midsurface.  The  r  is 
the  panel's  radius  of  curvature.  The  $'s  are  rotations  expressed 
as : 


* 


x 


<t>  =-w  +  (v/r)  (6) 

j  *y 

0  =0.5 (v  -u  ) 

» X*  tj 

Note  that  the  commas  in  the  subscripts  denote  partial 
differentiation  with  respect  to  the  parameter  that  follows  the 
comma • 

The  stresses  in  the  kth  lamina  can  be  expressed  in  terms  of 
the  midsurface  strains  and  curvatures  by  substituting  Eq.  (4) 
into  Eq.  (1)  (with  the  transformed  reduced  stiffness  matrix) 
resulting  in: 


’  ' 

a.. 
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rr  €° 

X 

K 

X 

Or 

J 
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The  resultant  forces  and  moments  per  unit  width  acting  on  a 


laminate  as  shown  in  Fig.  3  are  now  obtained  by  integrating  the 
stresses  in  each  lamina  through  the  laminate  thickness,  as: 


M  =  /  <TvZdZ 

x  J-t/2  x 

All  the  forces  and  moments  shown  in  Fig.  3  can  be  written  as 


where  N  is  the  total  number  of  laminae  and  and 

defined  in  Fig.  2. 


zk-l 


are 


One  can  now  substitute  Eq.  (7)  into  Eq.  (9)  to  arrive  at  an 
expression  for  the  force  and  moment  resultants  in  terms  of  the 
transformed  reduced  stiffness  matrix  and  the  midplane  strains  and 
curvatures.  Since  the  stiffness  matrix  is  constant  over  the 
lamina,  it  can  come  out  from  under  the  integration  sign.  Also 
note  that  the  midplane  strains  and  curvatures  are  independent  of 
Z;  therefore,  they  can  be  removed  from  under  the  integration  and 
summation  signs.  Thus,  the  force  and  moment  resultants  can  be 


written  as 


"11  "12  "16  • 
A21  A22  A26  ' 


B11  “12  “16  "x 

►  +  B21  B22  B2 6SA'J 


61  a62  a66  v  B6l  b62  b66  a  ■ 


B11  b12  Bl6lf€x  1  |D11  d12  Dl6lfAv 


=  B21  B22  B26  Irv  r  °21  °22  °26  N  As 
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where : 


Aij  <Zk  zk-i 


BiJ-°-5k?1(Qij,k  ,Zk_zk-i 


Dij  =0'3k?i  (Qj-j  >k 


The  A.  .  are  called  extensional  stiffnesses,  the  B.  .  are  called 
ij  ij 

coupling  stiffnesses,  and  the  D.  .  are  called  bending  stiffness. 

^  J 

It  is  noted  here  that  all  the  panels  considered  in  this 
thesis  have  symmetric  layups  about  the  midplane.  This  results  in 
no  coupling  between  extension  and  bending  since  the  B  matrix  is  a 


null  matrix  from  the  above  equations.  The  equations  for  A.  .  and 

**•  J 

are  used  in  subsequent  sections  to  explain  the  numerical 
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STAGSOl  Theory 


The  Structural  Analysis  of  General  Shells  (STAGSC-1) 
computer  code  [5,  6,  7]  is  a  finite  element  approach  based  on  the 
principle  of  minimum  potential  energy.  The  principle  of  minimum 
potential  energy  states  that  [11]  : 


"Among  all  admissible  configurations  of  a 
conservative  system,  those  that  satisfy  the  equations 
of  equilibrium  make  the  potential  energy  stationary 
with  respect  to  small  variations  of  displacement.  If 
the  stationary  condition  is  a  minimum,  the  equation 
state  is  stable." 


Loss  of  stability  (bifurcation)  occurs  when  the  second  variation 
of  a  systems’  total  potential  energy  ceases  to  be  positive 
definite.  The  following  discussion  outlines  how  STAGSC-1 
establishes  the  bifurcation  load  through  an  eigenvalue  technique. 
This  discussion  can  also  be  found  in  Ref.  [2], 

A  shell's  total  potential  energy  is  given  by 

V=U-W  (12) 

where  V  is  the  total  potential  energy,  U  is  the  internal  strain 


energy,  and  W  is  the  work  done  by  external  forces.  This  can  also 
be  written  as  [11] 

V.0.5  |d|T  [K1  |d|  -  jdjT  |r|  ,13, 
where  Jdj  is  nodal  degree  of  freedom  of  the  structure  vector,  [K] 
is  the  structure  stiffness  matrix,  |r|  is  the  structural  applied 


load  vector. 


For  a  composite  panel,  the  element  strain  energy  is  given  by 


[A]  [B] 

[B]  [D] 


Where  [A] ,  [B] ,  and  [D]  are  the  3x3  stiffness  matrices  given  by 
Eq.  (11),  the  j  t  j  is  the  midplane  strain  and  curvature  vector 
given  by  Eqs.  (5)  and  (6),  N  is  the  number  of  lamina  in  the 
laminate,  and  S  is  the  surface  area. 

Examination  of  Eqs.  (5)  and  (6)  reveals  that  the  strain 
vector  is  a  function  of  the  midplane  displacements,  their  first 


partial 

derivatives 

with  respect  to  x 

and 

Y> 

and  their  second 

partial 

derivatives 

of  w  with  respect 

to 

X 

and  y . 

Bauld  [9] 

carried 

out  the  integration  of  Eq. 

(14) 

and 

found 

that  this 

expression  for  the  strain  energy  consisted  of  three  distinct 
parts.  The  first  was  quadratic  in  displacements,  the  second  was 
cubic  in  displacements,  and  the  third  was  quartic  in 
displacements.  This  result  was  written  in  terms  of  the 
appropriate  element  shape  functions  and  nodal  degrees  of  freedom 
to  get  the  form  of  Eq.  (13)  .  The  result  for  the  total  potential 
energy  is  then  [9] 

V=(0*5Hqr+0,16N1qr+0-083N2qr)c^qr-Rrqr  (15) 

where  the  structure  nodal  degrees  of  freedom  are  replaced  by  the 
displacement  vector,  q.  is  the  system's  linear  stiffness 
matrix  with  no  dependence  on  the  displacement  vector,  q.  N1  ^  and 


14 


N2qr  are  matrices  with  linear  and  quadratic  dependence, 
respectively,  on  displacement. 

From  the  previous  statement  of  the  principle  of  total 
potential  energy,  setting  the  first  variation  of  Eq.  (15)  equal 
to  zero  gives: 

(Hqr+0 . 5Nlqr  +  0 . 3N2qr)qr-RQ  =0  ( 16 ) 

The  following  expression  is  found  by  satisfying  the 
condition  for  bifurcation  and  taking  the  second  variation  of  the 
total  potential  energy,  and  finding  the  point  where  it  is  no 
longer  positive  definate: 

DET  <Hqr+Nlqr+N2qr)=0  (17) 

Equation  (17)  is  used  by  STAGSC-1  to  solve  the  eigenvalue  problem 
of  the  form 

(H] +x [I] +\2 [J] =0  (18) 

where  [I]  and  [J]  represent  nonlinear  stiffness  matrices  in 
unknown  displacements  and  products  of  displacements, 
respectively.  For  a  linear  analysis , \  ,  is  the  proportionality 
constant  of  a  convenient  load  level  used  in  the  equilibrium  Eq. 
(16)  to  solve  for  the  unknown  displacements.  The  J  matrix,  which 
arises  from  the  prebuckling  rotations, is  often  omitted  as  it  will 
be  here.  The  quantities  H  and  I  or  equivalently,  H  and  N1  ^ 


are 


III.  Modelling 


Boundary  Conditions  and  Loading 

Figure  4  is  a  sketch  depicting  the  panel  analysed.  The 
boundary  conditions  were  such  that  a  simply  supported  panel  in 
pure  shear  was  modeled.  To  accomplish  this,  displacements  in  the 
v  direction  (circumferential)  were  allowed  to  move  on  all  sides 
)  except  the  bottom.  This  was  to  prevent  free  body  motion.  All 

other  displacements  remained  fixed.  The  only  rotations  allowed 
were  those  along  the  edges  in  the  direction  of  the  edge  and  all 
Z  rotations  normal  to  the  panel  (Rw) . 

The  loads  applied  as  shown  in  Fig.  4  were  line  loads.  That 
is,  they  have  units  of  force/unit  length.  As  such,  the  applied 
9  forces  NXy  and  N^x  were  equal  around  the  panel.  This  assured 

equilibrium  and  load  continuity  around  the  panel. 

With  this  combination  of  loads  and  boundary  conditions,  the 

9  only  load  that  was  actually  applied  to  the  panel  was  the  one 

along  the  top  of  the  panel.  This  is  due  to  the  fact  that  fixed 
boundry  conditions  override  applied  loads  in  the  STAGSC-1 
computer  code. 

Grid  Selection 

Both  references  [1]  and  [2]  indicate  that  for  accuracy  at 
least  five  nodes  should  be  present  in  each  half  sine  wave  of  the 
deformed  panel.  Initial  computer  runs  provided  data  indicating 
that  0.5  inch  grids  would  be  required  to  provide  these  five 
nodes.  This  grid  size  was  used  for  aspect  ratios  of  0.5  and  1.0 
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but  a  grid  size  of  1.0  inch  was  used  for  the  aspect  ratio  of  two. 
This  was  to  conserve  computer  time.  A  comparison  to  a  grid  size 
of  0.5  inch  proved  no  appreciable  loss  of  accuracy. 

Aspect  Ratios 

Three  aspect  ratios  were  investigated.  These  were  0.5,  1.0, 
and  2.0.  Aspect  ratio  is  defined  as  the  axial  length  (x)  divided 
by  the  circumferential  length  (y) .  The  circumferential  length 
remained  a  constant  12  inches  while  the  axial  length  was  6,  12, 
and  24  inches.  This  was  an  attempt  to  investigate  only  aspect 
ratio  effects  on  buckling  and  not  curvature  effects. 

Lamina  Orientations 

Three  lamina  orientations  were  investigated.  These  were 
[0,+45,-45,90]s  ,  [0,90J2s,  and  [+45,-45]2s.  These  provided  three 
symmetric  orientations  with  widely  varying  shear  coupling 
stiffness  terms  A-^ ,  A^ ,  ,  and  D2£  terms.  It  will  be  shown  in 
a  later  chapter  what  the  values  of  these  coefficients  are  and  how 
they  affect  the  buckling  load. 

Element  Selection 

STAGSC-1  has  two  elements  well  suited  for  modeling  a  curved 
panel,  the  QUAF  410  and  QUAF  411  elements  [6],  Though  Ref.  [2] 
indicates  the  QUAF  410  element  would  have  sufficed  for  this 
linear  study,  the  QUAF  411  element  was  selected  to  ensure 
accuracy.  This  element  and  the  degrees  of  freedom  are  shown  in 


The  additional  midside  tangential  displacement  degree  of 


freedom,  g..,  and  the  rotational  degree  of  freedom  of  the  edges 

iJ 

meeting  at  a  node, a  ,  make  the  411  element  better  suited  for 
nonlinear  analysis  of  cylindrical  panels  by  permitting  shear 
strain  matching  at  the  corners  [7] .  This  is  an  additional 
characteristic  the  410  element  does  not  posses  since  it  is  a  flat 
element.  Thus,  the  411  element  is  better  suited  to  a  cylindrical 
shell.  This  is  done  by  insuring  the  transverse  displacement  shape 
functions  are  of  the  same  order  as  the  in-plane  shape  functions. 
Therefore  the  in-plane  shape  functions  are  cubic  by  introducing 
the  first  normal  rotation, 0 „ ,  at  the  corners. 

Lt 


Material  Selection 

•  The  two  materials  used  in  this  analysis  are  graphite/epoxy 

and  kevlar/epoxy.  The  material  properties  for  the  graphite/epoxy 
were  the  same  used  in  Refs.  [2],  [3],  and  [4].  These  are: 

0  E1  =20500  ksi 

=  1300  ksi 


G12= 

750  ksi 

4 

k: 

i—1 

f\> 

ii 

0.33 

The 

material  properties  for 

kevlar/epoxy 

were 

found 

experimentally  in  Ref.  [14].  These  are: 

*  E-^  =10700  ksi 

E2  =  600  ksi 

G12=  ksi 

*  *  i °.25 


% 
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Flat  Plate  Check 


To  check  to  see  if  the  boundary  conditions  and  loading 
geometry  represent  a  simply  suported  panel  in  pure  shear,  a 
comparision  was  carried  out  of  the  results  arrived  at  with 
STAGSC-1  with  the  analytic  results  of  a  flat  square  panel  [12]. 
The  analytic  solution  of  the  flat  plate  buckling  is  a  Galerkin 
series  solution. 

To  match  the  stiffness  criterion  given  for  the  analytic 
solution,  the  material  properties  had  to  be  changed  to: 

E^  =20500  ksi 
E2  =  2050  ksi 

G12=  ^70  ksi 
V- ~=  0.33 

Also,  all  eight  plies  of  the  panel  were  changed  to  a  zero 
orientation  since  the  analysis  was  on  an  orthotropic  plate. 
Finally,  the  radius  of  curvature  of  the  panel  was  increased  to  24 
inches  to  better  simulate  a  flat  plate.  All  other  parameters  were 
kept  the  same  as  in  the  baseline  panel  with  aspect  ratio  1.0. 

The  Galerkin  method  arrived  at  a  buckling  load  for  the  flat 
plate  of  NXy=23.7  #/in.  The  STAGSC-1  analysis  obtained  a  buckling 
value  of  NXy=60  #/in.  The  diferrence  in  the  predicted  loads  can 
be  attributed  to  two  things.  First,  even  though  the  radius  of 
curvature  was  increased  to  24  inches,  the  curved  panel  should  be 
able  to  carry  more  load  then  the  flat  plate.  Secondly,  the  flat 
plate  analysis  does  not  take  into  consideration  the  membrane 
stresses  that  the  STAGSC-1  analysis  does.  Under  the  applied  shear 
load,  these  membrane  stresses  will  be  tension  thus  producing  a 


22 


higher  load  carrying  capability. 

It  is  shown  in  chapter  4  that  a  curved  composite  panel  does 
not  react  as  a  flat  plate  does  under  pure  shear.  Looking  back  at 
this  analysis  shows  that  24  inches  is  sufficient  curvature  to 
prevent  a  pure  shear  reaction.  Hence,  these  results  are 
misleading.  The  next  check  better  approximates  a  flat  plate  with 
a  larger  radius  of  curvature. 


Curved  Panel  Check 

Similar  to  the  flat  plate  check,  a  comparison  of  the 
STAGSC-1  analysis  was  made  to  a  curved  panel  analytical  solution 
to  check  boundary  conditions  and  loading  geometry  [13].  This  was 
also  a  Galerkin  method  solution. 

The  material  properties  had  to  again  be  modified  to  meet  the 
analytic  analysis.  This  time  the  material  properties  are: 

Ex  =20500  ksi 
E2  =  1500  ksi 
G^2=  730  ksi 

V12=  0,25 

The  radius  of  curvature  of  the  panel  also  had  to  be  increased  to 
40  inches  to  avoid  extrapolating  the  documented  results.  All 
other  parameters  were  the  same  as  the  [+45,-45]2s  baseline  panel 
with  an  aspect  ratio  of  one. 

The  analytical  Galerkin  method  gave  a  bifurcation  applied 
stress  of  a  =3000  psi.  The  STAGSC-1  analysis  again  gave  a  higher 
value  of  <T  =3650  psi.  This  difference  can  again  be  explained  by 
the  Galerkin  method  not  taking  membrane  stresses  into  account. 


The  bifurcation  eigenvectors  and  stress  distribution  for 
this  panel  are  shown  in  Fig.  6.  The  eigenvectors  are  oriented  at 
45°  as  would  be  expected  for  a  flat  plate  loaded  under  pure 
shear.  The  stress  distribution  also  reflects  a  panel  in  pure 
shear.  This  indicates  the  boundary  and  loading  conditions  are 
adequate  to  model  a  flat  plate  in  pure  shear.  That  is,  a  curved 
panel  with  a  radius  of  curvature  of  40  inches  can  be  considered  a 
flat  plate. 


Figure  6.  Eigenvectors  for  and  .Stress  distribution 

Throughout  §+5]  2s  l  anel  with  40  Inch  Radius 
of  Curvature.  Cots  Represent  Negative 
Radial  Displacements 


IV.  Results 


Numerical  Analysis  Results 

The  ASD  Cyber  computer  system  was  used  to  implement  the 
STAGSC-1  computer  code.  A  sample  of  the  input  decks  used  is  shown 
in  Appendix  A.  Table  1  shows  the  bifurcation  load  computed  for 
all  the  panel  configurations  studied. 

The  eigenvalue  routine  that  the  STAGSC-1  code  uses  computes 
all  eigenvalues,  both  positve  and  negative.  The  values  in  Table  1 
are  the  smallest  positive  values.  These  correspond  to  the  first 
bifurcation  load.  The  negative  values  have  no  physical 
significance.  They  are  merely  mathematical  anomales.  It  was  noted 
in  the  computer  output  that  a  negative  eigenvalue  normally 
occurred  with  approximately  the  same  magnitude  as  each  positive 
eigenvalue . 

Examination  of  Table  1  shows  a  few  general  trends.  First, 

the  results  for  graphite-epoxy  are  consistantly  higher  than  those 

for  the  Kevlar-epoxy .  This  is  directly  attributable  to  the  higher 

material  stiffness  properties  of  the  graphite-epoxy  over  the 

Kevlar-epoxy.  Secondly,  the  bifurcation  load  decreases  with 

increasing  aspect  ratio.  This  is  due  to  the  increase  in  panel 

area  with  the  increasing  area  ratio.  It  is  physically  obvious 

that  an  infinitly  larger  panel  will  not  be  able  to  carry  a  large 

line  load,  N,  ,  since  the  side  length  qoes  to  infinity, 
xy 

Similarlly,  an  infinitly  small  panel  will  be  able  to  carry  a  very 
large  line  load,  ,  since  the  side  length  goes  to  zero.  This  is 
in  line  with  the  results  found  in  Ref.  [12].  Lastly,  the 
bifurcation  load  for  the  [-*-45,— 451^3  is  the  highest  with  the 
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[0 , +45 ,-45 , 9  0 ] s  next  and  the  [0,90]^  panel  being  the  weakest. 
This  corresponds  to  the  order  of  the  largest  to  smallest  value 
for  the  and  stiffness  terms  (see  Appendix  B) .  This  is 
expected  due  to  the  loading  applied  to  the  panels. 

The  eigenvectors  are  represented  graphically  in  Figs.  7  thru 
12.  These  figures  are  2-dimensional  representations  of  the  curved 
panels  in  which  the  panels  have  been  flattened  out.  The  dots  in 
these  figures  represent  negative  radial  displacements.  These 
figures  clearly  show  the  troughs  in  the  bifurcated  panel. 


Eigenvectors 

The  eigenvectors  shown  in  Figs.  7  thru  12  do  not  display  the 
expected  results  of  a  panel  loaded  in  pure  shear.  The  expected 
results  would  have  the  displacement  troughs  running  diagonally 
across  the  entire  element  like  those  shown  in  Fig  6.  The  computed 
eigenvectors  show  this  trend  for  a  portion  of  the  panel  but  not 
for  the  entire  height.  This  is  attributed  to  the  shell  analysis 
of  the  curved  panels  resulting  in  different  internal  stress 
distributions  than  those  for  a  flat  plate. 

This  different  stress  distribution  is  shown  in  Fig.  13. 
These  zones  of  resultant  stress  distributions  were  found  to  be 


consistant  in  all  the  panel  configurations  examined.  Figure  14 
shows  the  principal  stress  conversions  from  the  various  computed 
stress  states  showen  in  Fig.  13.  These  principal  stresses  explain 
most  of  the  eigenvectors  showen  in  Figs.  7  thru  12.  This  is  done 
keeping  in  mind  the  troughs  will  run  perpendicular  to  the 
compressive  load  line  of  action.  Two  zones,  however,  are  not 
explained  by  this. 


Figure  7*  Eigenvectors  for  Graphite/Epoxy  [0,t45,90]  •  Dots 

Represent  Negative  Radial  Displacements. 
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or 


Aspect  Ratio-1.0 


Aspect  Ratio-2.0 


Figure  10.  Eigenvectors  for  Kevlar/h'poxy  £o,t 45 ,90)  . 

Represent  Negative  Radial  Displacements . 


aspect  Ratio-1.0  aspect  Ratio-2.0 


Figure  13*  Nx»  N^.,  and  N  Distribution  Throughout  iranels 
from  Computer  Analysis 


These  two  zones,  the  lower  half  of  both  the  panels  of  aspect 
ratio  one  and  two,  can  be  explained  by  referring  to  Fig.  15.  This 
figure  shows  the  prebifurcation  radial  displacements  of  the 


graphite/epoxy  [0 , +45 ,-45 , 9 0 ] s  panel  for  an  example.  Since  Fig. 
13  shows  no  stresses  in  the  lower  portion  of  the  panel  with  an 
aspect  ratio  of  two,  the  radial  displacements  after  bifurcation 
will  be  the  same  as  the  radial  displacement  before  bifurcation. 
This  is  the  case  on  comparison  of  Figs.  7  and  15. 

Explaining  the  lower  portion  of  the  panels  with  an  aspect 
ratio  of  one  is  not  as  clear.  Figures  13  and  14  would  indicate 
the  eigenvector  should  be  oriented  at  some  angle  to  the 
horizontal  other  than  the  90°  indicated  in  Fig.  7.  This  must  be 
due  to  an  interaction  between  the  prebifurcation  radial 
displacements  and  those  that  occur  after  bifurcation.  It  was 
observed  that  the  prebifurcation  radial  displacements  were,  in 
general,  small  in  comparison  to  the  shape  of  the  bifurcation 
displacements.  This  generalized  comparison  is  not  true  for  the 
lower  portion  of  the  panels  of  aspect  ratio  one.  This  is 
attributed  to  the  decreasing  magnitude  of  NXy  in  this  region. 
Therefore,  the  bifurcation  eigenvector  will  be  made  up  of  a 
combination  of  the  prebuckeled  and  bifurcation  radial 
displacements . 

The  different  stress  distributions  for  the  different  aspect 
ratios  can  be  qualitativly  explained  by  considering  the  boundary 
conditions  used  and  the  applied  load.  Since  the  only  load  that  is 
actually  applied  to  the  panel  is  the  line  load  along  the  top  of 
the  curved  panel,  a  twist  is  applied  to  the  panel.  This  twist 


Aspect  Ratio-1.0 


Aspect  Ratio-2.0 


Figure  15-  1  rebifurcation  Radial  Displacements  of  Graphite/Dpoxy 
[0,^45,90]  _  lanel.  Dots  Represent  Negative 

“  5 

Radial  Displacements. 
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results  in  the  panel  trying  to  warp.  This  is  something  that  is 
not  seen  with  flat  panels.  This  warping  is  the  panel  trying  to 
equilibriate  the  N  resultant  force.  Since  the  boundary 

A 

conditions  in  this  analysis  restricted  vertical  displacements  on 
the  vertical  edges,  this  warping  was  restricted  and  increasses 
the  membrane  stresses.  This  results  in  the  resultant  N  force 

A 

throughout  the  panels. 

The  reason  the  aspect  ratio  effects  this  N  distribution  is 

A 

the  coupling  between  the  vertical  side  boundary  conditions  and 
the  applied  load  and  the  fixed  bottom.  The  panel  with  an  aspect 
ratio  of  one  has  the  greatest  distribution  of  the  Nx  resultant 
force.  This  is  due  to  the  vertical  boundaries  and  the  top 
load/fixed  bottom  being  in  closest  proximity.  As  the  aspect  ratio 
decreasses  to  0.5,  the  top  load  and  bottom  move  closer  together 
but  the  vertical  boundaries  move  proportionatly  further  appart. 
The  converse  is  true  when  the  aspect  ratio  increases  to  two.  This 
indicates  a  smaller  resultant  N  force  influence  on  the  aspect 

A 

ratio  0.5  and  2.0  then  with  an  aspect  ratio  of  1.0  which  is  the 
case . 

The  above  discussion  verifies  that  a  curved  panel  analysed 
as  a  shell  and  loaded  in  pure  shear  do  not  react  as  a  flat  plate 
subjected  to  the  same  load  arrangement. 

Nondimensionalization  Analysis 

Material  Analysis.  The  first  analysis  was  done  to  remove 
the  material  property  dependence  on  the  bifurcation  load.  It  was 
expected  that  the  shear  modulus  of  elasticity,  , 


would  be 


useful  in  accomplishing  this.  As  such,  the  bifurcation  load 


(remember,  this  load  is  dimensionally,  force/unit  length)  was 
divided  by  G-^  and  the  circumferential  length  of  12  inches.  The 
circumferential  length  was  used  instead  of  the  variable  height 
because  the  resulting  value  is  plotted  against  the  aspect  ratio. 
These  results  are  shown  in  Fig.  16. 

The  resulting  plots  show  clearly  that  this  assumption 


relative  to  G-^  is  in  error.  The  curves  for  the  two  materials  do 
not  display  sufficient  correlation  to  justify  the  conclusion  that 
the  shear  modulus,  is  the  proper  material  property  for  this 
analysis . 

Instead  of  using  G^.,,  the  extensional  modulus  of  elasticity, 
E^  was  used.  Again,  the  bifurcation  load  was  divided  by  and 
the  circumferential  length.  The  results  of  this  analysis  are 
shown  in  Fig.  17.  This  analysis  provided  a  relationship  with 
reasonably  close  correlation. 

The  only  area  where  this  correlation  is  not  very  close  is 
with  the  [0,90]2S  panel  with  aspect  ratio  greater  than  one.  The 
load  carrying  capability  of  this  panel  with  an  aspect  ratio  of 
two  is  about  the  same  as  that  of  a  panel  with  an  aspect  ratio  of 
one.  This  can  be  explained  by  referring  back  to  Fig.  13  and  the 
stiffness  coief f icints ,  and  <3^ven  in  Appendix  B.  Since 
Dl6  anc*  °2 6  have  values  of  zero  for  this  lamina  orientation, 
these  panels  can  only  resist  shear  along  a  45°  diagonal.  However, 


these  panels  have  great  stiffness  in  the  axial  directions.  Thus, 
the  bifurcation  of  these  panels  will  be  due  to  the  shear  load  and 
not  the  axial  components  of  load.  Keeping  this  in  mind,  Fig.  13 
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shows  a  very  similar  relationship  of  stress  distribution  over  the 
same  total  area  between  the  panels  with  an  aspect  ratio  of  one 
and  two.  Hence,  the  bifurcation  load  should  be  similar. 

It  is  also  apparent  in  Fig.  17  that  the  curves  are  not 
linear.  This  is  again  due  to  the  fact  that  as  aspect  ratio 
decreases,  the  line  load  carrying  capability  will  go  to  infinity 
since  edge  length  is  approaching  zero.  Converslly,  as  the  aspect 
ratio  increases,  edge  length  goes  to  infinity  and  thus  the  line 
load  carrying  capabiliy  goes  to  zero.  Furthermore,  the  various 
membrane  resultant  forces  are  affecting  the  linear  relations  of 
load  versus  aspect  ratios. 

This  technique  of  using  the  Young's  modulus  to  remove 
material  property  dependence  from  buckling  analysis  is  also  used 
in  Ref.  [13].  In  this  reference  the  transverse  modulus,  Eo  ,  was 
used  as  opposed  to  the  longitudnal  modulus,  E-^  ,  as  used  here. 
Examining  the  moduli  of  the  two  materials  used  in  this  thesis 
reveals  that  E-^  (Gr-Ep)  /E-^  (Ke-Ep)  and  E2  (Gr-Ep)  (Gr-Ep)  are  the 
same.  Hence,  Fig.  17  would  have  the  same  trend,  just  different 
magnitude. 

Qualitativly ,  there  is  a  reason  for  the  extensional  modulus 
being  the  important  maerial  property  for  shear  buckling  analysis 
for  these  panels  when  shear  buckling  is  considered  in  isotropic 
materials.  As  shown  in  Fig.  14,  the  failure  mode  for  shear 
loading  is  actually  a  compresive  failure  i  an  off  axis  direction. 
For  an  isotropic  material  this  compresion  faiure  is  resisted  by 
E.  Since  the  [ 0  , + 4 5 , -45 , 90 ] 0  panel  is  quasi-isotropic,  it  will 

O 

also  have  this  trait.  Looking  back  at  Fig.  17  showes  that  this 
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panel  results  in  the  best  correlation.  The  other  two  panels 
stiffness  properties  move  further  away  from  quasi-isotropic  and 
thus  have  a  lower  degree  of  correation. 

Aspect  Ratio  Analysis .  This  analysis  was  accomplished  by 
normalizing  the  bifurcation  loads  to  the  bifurcation  loads  of  the 
panels  with  an  aspect  ratio  of  one.  In  Fig.  18  these  normalized 
loads  are  plotted  against  aspect  ratios.  The  result  is  a 
relationship  having  all  the  lamina  orientations  in  close 
agreement  with  each  other,  even  between  materials.  The  only  point 
out  of  agreement  is  the  [0,903^2  panel  above  an  aspect  ratio  of 
one.  This  can  again  be  attributed  to,  as  in  the  discussion  on 
material  properties,  the  fact  that  the  failure  mode  is  the  same 
over  the  same  total  area  for  panels  of  both  an  aspect  ratio  of 
one  and  two.  This  again  results  in  a  bifurcation  load  for  these 
two  aspect  ratios  to  be  very  similar. 

Also,  the  curves  in  Fig.  18  are  nonlinear  for  the  same 
reason  as  given  in  the  discussion  on  material  properties  and  will 


not  be  repeated  here. 

An  explanation  of  the  close  correlation  shown  in  Fig.  18  is 
that  the  percentage  change  in  bifurcation  load  in  any  ply  lay-up 
is  constant  no  matter  what  ply  layup  is  being  considered.  This  is 
obviously  not  the  case  with  the  t0,90]2s  orientation  for  the 
reasons  outlined  above.  A  more  complete  understanding  will 
require  more  lay-ups  to  be  analysed. 

Lamina  Orientation  Analysis.  The  last  parameter  to  be 
studied  is  the  effect  of  lamina  orientation  on  bifurcation  load. 
This  analysis  again  normalized  all  the  bifurcation  loads  to  a 


standard  load.  This  time  the  bifurcation  loads  were  normalized  to 


the  bifurcation  loads  of  the  [0,90] panels.  The  parameter  these 
normalized  loads  were  plotted  against  were  the  coupling  bending 
stiffness  terms  and  £>26' 

These  stiffness  terms  were  used  because  they  are  the  ones 
that  mostly  resist  shear  buckling  loads  [10]  .  The  values  for 
these  stiffness  terms  can  be  found  in  Appendix  B  for  all  the 
lamina  orientations. 

The  results  of  this  analysis  are  shown  in  Fig.  19.  The 
trends  are  linear  with  the  dependence  on  aspect  ratios  0.5  and 
1.0  removed.  The  panels  with  an  aspect  ratio  of  two  do  not  follow 
this  trend.  This  is  due  again  to  the  fact  that  these  panels  do 
not  carry  all  their  load  in  pure  shear.  Since  this  is  the  case, 
the  bending  stiffness  terms,  and  D ^  are  not  the  values  these 
bifurcation  loads  should  be  plotted  against.  Instead,  a 
combination  of  bending  stiffnesses  ,  and  as  well  as 
axial  extension  A-,  ^  ,  A„x  ,  and  A^ ✓  should  be  used. 
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Figure  19*  Effect  of  Lamina  Orientation  on  Bifurcation  Load 
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V.  Conclusions 

Based  on  the  analysis  completed  in  this  thesis  the  following 
conclusions  have  been  made: 

1.  The  boundary  conditions  used  in  this  thesis  are  valid 
for  applying  a  pure  shear  load  to  a  flat  plate  for  a  STAGSC-1 
computer  analysis. 

2.  Curved  panels  that  are  analysed  as  shells  loaded  under 
pure  shear  develop  an  interior  shear  distribution  unlike  that  of 
flat  plates  under  the  same  load  configuration. 

3.  The  Young's  modulus,  E-^  or  E2  ,  is  a  better  indicator  of 
shear  stiffness  than  the  shear  modulus,  G-^  for  the  composite 
panels  analysed. 

4.  As  aspect  ratios  approach  zero,  the  bifurcation  load 
approaches  infinity.  Converslly,  as  the  aspect  ratio  approaces 
infinity,  the  bifurcation  load  generally  approaches  zero. 

5.  A  panel  with  a  ply  orientation  of  [0,90] 2s  has 
approximately  a  constant  load  carrying  capability  above  an  aspect 
ratio  of  one. 

6.  A  panel's  aspect  ratio  effects  the  interior  distribution 
of  stresses. 

7.  The  bending  stiffness  coefficients,  and  D2£  ,  effect 
the  load  carrying  capability  of  a  curved  composite  panel  as  they 
do  for  flat  plates.  This  effect  results  in  a  higher  load  carrying 
capability  for  higher  values  of  and  D2£  . 
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Appendix  A 


Sample  Input 

An  example  of  an  input  deck  used  in  this  thesis  is  shown  in 
Fig.  1-A.  A  detained  discription  of  each  input  line  is  given  in 
Ref.  6.  Special  caution  must  be  taken  when  inputing  Poission's 
ratio.  Most  material  property  tables  give 
code  uses  . 


but  the  STAGSC-1 


LINEAR  SHEAR  ANALYSIS,  ASPECT  RATIO  1.0,  GR-EP,  0, +45, -45, 90 

1  $  B-l  LINEAR  STATIC  ANALYSIS 

10$  B-2  1  SHELL  0  ELEMENTS 

1010$  B-3  1  MATERIAL  1  SHELL  WALL 

1.0  $  C-l  LOAD  MULTIPLIER  RECORD 

1  0  650  $  D-2  EIGENVALUE  CONTROL  CARD 

5  $  D-3  CLUSTER  DEFINITION  RECORD 

25  25  $  F-l  25  ROWS  25  COLUMNS 

1  0  $  1-1  MATERIAL  NUMBER 

20.5E06  0.0212429  .75E06  0.0  1.0  1.3E06  1.0  $  1-2  MATERIAL 
PROPERTIES 

118$  K-l  WALL  CONFIGURATION 

1  .005  0.0  $  K-2  PLY  MATERIAL,  THICKNESS,  ORIENTATION 
1  .005  45.0  $  K-2 
1  .005  -45.0  $  K-2 
1  .005  90.0  $  K-2 
1  .005  90.0  $  K-2 
1  .005  -45.0  $  K-2 
1  .005  45.0  $  K-2 

1  .005  0.0  $  K-2  END  OF  ALL  8  PLYS 

5  $  M— 1  CYLINDRICAL  SHELL  UNIT 

0.0  12.0  0.0  57.296  12.0  $  M-2A  SHELL  GEOMETRY 

1  0  $  M-5  SHELL  WALL  RECORD 

411  $  N-l  ELEMENT  NUMBER 

0000$  P-1  BOUNDRY  CONDITIONS  SPECIFIED  ON  P-2  RECORD 

010  011  $  P- :  V,  RV,  RW  FREE  ON  TOP 

010  101  $  P-2  V,  RU,  RW  FREE  RT  SIDE 

000  011  $  P-2  RV,  BW  FREE  ON  BOTTOM 

010  101  $  P-2  V,  RU  , RW  FREE  LT  SIDE 

1  $  Q-l  1  LOAD  SYSTEM 

140$  Q-2  1  LOAD  SET  SYSTEM  A 

1.0  2  2  1  $  Q-3  TOP  LINE  LOAD 

-1.0  3  0  0  13  $  Q-3  RIGHT  LINE  LOAD 

-1.0  2  2  13  $  Q-3  BOTTOM  LINE  LOAD 

1.0  3  1  0  1  $  Q-3  LEFT  LINE  LOAD 

1100$  R-l  PRINT  OUTPUT 


Figure  1-A.  Sample  STAGSC-1  Input 


Appendix  B 


Stiffness  Coeificients 

The  following  are  the  extensianal  and  bending  stiffness 
coeifficient  matrices  for  the  panels  investigated  in  this  thesis. 
The  equations  given  in  chapter  two  were  used  with  the  material 
properties  given  in  chapter  three.  Note  that  since  the  panels  are 
all  symmetric,  the  coupling  stiffness  matrix  is  the  null  matrix. 

Graphite-Epoxy 
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A= 


0.35  0.11  0.00 

77.0  12.0  4.80 

0.11  0.35  0.00 
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0.00  0.00  4.00 
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[+45 ,-45] 2S 


■0.26  0.20  0.00' 

T4.0  26.0  9.61T 

A= 

0.20  0.26  0.00 

X  106  #/  in  D= 

26.0  34.0  9.60 

0.00  0.00  0.21 

9.60  9.60  28.0 

Kevlar-Epoxy 
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